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Abstract 

We perform for the Large Hadron Collider (LHC) a detailed study of charged Higgs 
boson production via the top-bottom quark associated mode followed by decays into 
a chargino and a neutralino, with masses and couplings as given by the general Min- 
imal Supersymmetric Standard Model (MSSM). We focus our attention on the region 
of parameter space with mjj± > mt and intermediate values of tan/3, where identifica- 
tion of via decays into Standard Model (SM) particles has proven to be ineffective. 
Modelling the CMS detector, we find that a signature consisting of three hard leptons 
accompanied by a hadronically reconstructed top quark plus substantial missing trans- 
verse energy, which may result from — > Xi^2^i 234 decays, can be made viable over 
a large variety of initially overwhelming SM and MSSM backgrounds, provided MSSM 
input parameters are favourable: notably, small |^| and light sleptons are important 
prerequisites. We quantify these statements by performing a fairly extensive scan of 
the parameter space, including realistic hadron-level simulations, and delineate some 
potential discovery regions. 
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1 Introduction 



A pair of spin-less charged Higgs bosons, (with mass m_H-±), arises in any Two-Higgs 
Doublet Model (2HDM) alongside a trio of neutral Higgs bosons — the CP-even 'light' h 
and 'heavy' H {i.e., with rrih < tuh) scalars and the CP-odd pseudoscalar ^4 (with mass 
itia)- Embedding a Type II 2HDM inside the attractive theoretical framework provided by 
Supersymmetry (SUSY) yields the MSSM (see [1]), wherein the particle content is hmited to 
the known SM states (fermions and gauge bosons), their 'sparticle' counterparts (sfermions and 
gauginos) plus the five aforementioned Higgs bosons and their respective Higgsinos. Among 
the new massive sparticles predicted in the MSSM are the charginos and the neutralinos^, 
which are the mass eigenstate mixtures of the electroweak (EW) gauginos and the Higgsinos. 
Previous papers [2, 3] have demonstrated that decays into a chargino and a neutralino 
can probe regions of the MSSM parameter space where charged Higgs boson decays into SM 
particles and other Higgs bosons are swamped by backgrounds. In particular, tan (3 (the ratio 
of the vacuum expectation values of the up- type and down-type Higgs doublets) values between 
3 and 10 were found to be in part accessible via 

X1X23 decay modes {i.e., charged 
Higgs boson decays into the hghtest chargino and the second or third heaviest neutrahno), 
when the final state includes three leptons (meaning electrons and/or muons)^. 

Such tan (3 values fall in the so-called 'intermediate' regime wherein decays to SM 
objects (which may include neutral MSSM Higgs bosons) are undetectable at the LHC irre- 
spective of the values chosen for other MSSM input parameters^. This zone of undetectability, 
in part due to the ~ {ml tan P'^+ml/ tan /3^) couphng of the main pp — > tH~X + c.c. produc- 
tion mode"^, begins around tan/3 = 6 or 7 for m/f± ~ rrit and spreads to encompass more and 
more tan f3 values (say, between 3 and 20) as mH± grows larger. The rate suppression may be 
further exacerbated by the same tan d dependence in the H~ hi decays if there are other 
competing decay channels — naturally, if the — > ht branching ratio (BR) is ~ 1, it will 
remain so and there is no additional suppression because of the bottom-top decay rate. The 
alternative MSSM decay channel H~ hW~ , which also yields bbW~ intermediate states 
(since h ^ bb), is only relevant within a minuscule tan/3 interval (roughly tan/3 ^2 — 3) for 
mH± i,mt — this lies close to the LEP2's excluded region. Then there is H~ tVt, which is 
limited to larger tan/? values^, at best offering coverage down to tan/3 ~ 10 for mj^± ~ rrit 
and contracting to even higher tan/3 values as mH± grows larger [7]. (See references in [2, 3] 

^We will refer to the charginos and neutralinos collectively as 'inos'. 

^Thc process is further identified by a hadronically reconstructed top quark from the tH^X (or tH^X) 
production process, and via substantial missing transverse momentum from the lightest neutralinos, XiS, the 
stable Lightest Supersymmetric Particles (LSPs) which must eventually result from decays of the inos. 

■'Not coincidcntally, in roughly the same area coverage via the neutral Higgs sector is questionable [4, 5], 
particularly if the integrated luminosity is limited (say, ~30fb^^). In fact, the one neutral Higgs boson that 
may be detectable typically mimics a SM Higgs boson (this is the so-called 'decoupling scenario'). 

"^Charged Higgs bosons might also be produced in sparticle decays [6]. In particular, decays of gluinos 
and/or squarks could lead to copious numbers of if containing events if these strongly-interacting sparticles 
are light enough to be abundantly produced. However, such events may also fail to contain a top quark or 
possess an excessive number of jets and so not satisfy our signal requirements. Here we neglect such production 
processes. 

^The H~ — > sc mode has a much reduced scope in comparison, because of the large QCD background. 
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for a list of phenomenological analyses of these SM decay modes of a charged Higgs boson.) 

Considering such limitations, it is worthwhile pursuing further the i?^ — > inos decay modes 
initially probed in [2, 3]^, expanding upon the results found therein and placing the analysis 
in a sounder phenomenological context. The improvements found herein go in three general 
directions. Firstly, the allowable parameter space is covered far more thoroughly, incorporating 
all possible chargino-neutralino decay modes into the analysis and including every conceivable 
path leading from a charged Higgs boson to a three leptons plus invisible energy final state. 
Secondly, investigation of the role of on- and off-shell sleptons (the SUSY partners of the 
leptons) is considerably deepened: as noted in the previous studies, if there is a light slepton, 
the leptonic BRs of the inos can be significantly enhanced (especially those of X2 and/or Xs)- 
Thirdly, signals are herein studied within a full event generator environment modelling the 
CMS detector and also includes an improved background analysis that encompasses potential 
MSSM background processes ([3] was a very prehminary account in both these respects while 
[2] only considered SM backgrounds and was carried out solely at the parton level). 

The legacy of the CERN e+e~ collider is a model independent limit on from charged 
Higgs pair production of order Mw± — 78.6 GeV is the current LEP2 bound [11]. Further, 
the current lower Higgs boson mass bound of approximately 114 GeV [11] can be converted 
within the MSSM into a minimal value for mH± of ~130-140GeV, for tan/3 ~ 3-4. This 
bound grows rapidly stronger as tan/9 is decreased while tapering very gradually as tan/9 is 
increased (staying in the 110-125 GeV interval for tan/9 > 6). For m^^i < rrit. charged Higgs 
bosons could be discovered during Run 2 of the FNAL Tevatron [12], which has ah^cady begun 
taking data at ^/s^p — 2 TeV, by exploiting their production in top and antitop quark decays 
{t — > bH'^ + c.c.) followed by H~ t~Vt -\- c.c. [13]. In contrast, if mH± S rrit (our definition 
of a 'heavy' charged Higgs boson), one will necessarily have to wait until the advent of the 
LHC at CERN, with \/Spp = 14 TeV, and thus this study will concentrate upon charged Higgs 
boson masses well above that of the top (anti) quark. This will also provide ample phase space 
to allow for decays into sparticles with masses above current experimental bounds. 

There are also other processes where charged Higgs bosons (or A, to whose mass that of 
the is closely tied) enter as virtual particles at the one-loop level. These include neutral 
meson mixing {K^K^, D°D° or B^B°) and Z ^ hi (Rb) [14], 6 ^ 57 decays [14, 15], b crUr 
decays [16] and the anomalous muon magnetic dipole moment [17]. The b ^ sj decays are 
generally thought to be the most constraining [14] (6 — > crPr becomes significant for very 
high values of tan/3). In the MSSM, one-loop diagrams with either a H^-t loop or a Xi~i 
can give meaningful contributions to 6 — > 57 decay processes (loops involving gluinos are not 
significant if gluinos are relatively heavy), and these two contributions may come either with 
the same or with opposite signs. Thus 6 — > 57 restrictions on uih^ are linked to the masses of 
the lighter chargino and the stops. The b ^ sj decays and the other higher order processes 
may well exclude some regions of the MSSM parameter space that are still allowed by the more 
direct limits from Higgs boson and sparticle searches at LEP2. However, definite bounds are 
quite difficult to delineate without restricting oneself to some subset of the allowed parameter 
space of the general MSSM by specifying a mechanism for how SUSY is to be broken (and 

^Hadron collider signals from neutral MSSM Higgs boson decays into inos were studied in [5, 8], while 
MSSM Higgs bosons BRs to inos, emphasising invisible decays to a pair of LSPs, were presented in [9, 10]. 



2 



in general linking together what in the general MSSM are independent input parameters). 
Studies which have delineated excluded regions resulting from these processes have invariably 
included additional assumptions about the behaviour of the theory at higher energy scales 
— such as in Minimal Supergravity (mSUGRA) for example, for which next-to-leading order 
(NLO) calculations have recently been performed [15]. There are also significant uncertainties 
in translating the experimental results into clear predictions about MSSM parameters [18]. 
Concerning limits from recent {g — 2)^ measurements, these are most restrictive [17] when 
tan (3 is low (< 3) - a case which is not of particular interest for our process - and may be 
relaxed when smuons are light - a case which is of particular interest for our process. 

2 MSSM Parameter Space 

Analysing the usefulness of chargino-neutralino decays within the general MSSM pa- 

rameter space is a fairly involved undertaking since many independent input parameters asso- 
ciated with just about all the (s)particle sectors of the model can play crucial roles. Prom the 
Higgs sector we of course have tan /3 along with one input Higgs boson mass, taken herein as 
rriA, to which the tree-level masses of all the other Higgs bosons are pegged. These two inputs 
are largely sufficient for the SM decay modes, assuming sparticlc decay modes arc closed. 

Squark masses, particularly stop masses, can drive significant radiative corrections to the 
tree-level Higgs boson masses, especially to m^. In contrast, higher order corrections to the 
tree-level relation = m\ + M^± are typically quite small [19]. Thus the signal rate is 
insensitive to the choice of squark-sector inputs. Nevertheless, the coloured-sparticle sector 
affects the analysis in peripheral — but potentially crucial — ways. Firstly, the choice of 
the stop mass inputs can affect what regions of the MSSM parameter space are excluded via 
Higgstrahlung or the aforementioned 6 — > 57 processes. The former would suggest choosing 
high stop inputs to help push rrih up above the LEP2 bounds, while the latter might prefer 
low stop inputs to cancel corrections due to a light chargino. Be such arguments as they 
may, there is considerable uncertainty in the resulting limits on the general MSSM parameter 
space, and these issues will not be addressed further. The second consideration is the size of 
squark and gluino backgrounds to our signature. Discussion of this will be postponed until 
the end of this section. 

To specify the ino sector, the parameters M2 and /i, in addition to tan f3, are required. Mi 
is assumed to be determined from M2 via gaugino unification (i.e., Mi = ^tan^ 6wM2). This 
will determine the tree-level masses (to which the radiative corrections are quite modest) of the 
inos along with their couplings to the Higgs bosons. However, this is not enough, for the inos 
(except for Xi) must also decay - preferably into leptons for easy detection. To calculate the 
leptonic ino BRs. one must designate the properties of the slepton sector, since light sleptons 
can greatly enhance said BRs [2, 3, 20]. Inputs (assumed to be fiavour-diagonal) from the 
slepton sector are the left and right soft slepton masses for each of the three generations 
(selectrons, smuons, and staus) and the trilinear '^-terms' which come attached to Yukawa 
factors and thus only A^- has a potential impact. A priori, all six left and right mass inputs 
(and At) are independent. However, in most models currently advocated, one has mg^ ~ m^ji^ 
and rrie^ ~ m^^. We will assume such equalities to hold. 
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Figure 1: Minimum allowed soft slepton mass, given constraints as described in the text. In 
the upper (lower) two plots, soft stau mass inputs are set 100 GeV above (degenerate with) 
those of the first two generations. = in all cases. The shaded areas are excluded by LEP. 



To maximise leptonic ino BR enhancement, sleptons should be made as light as possible. 
But direct searches at LEP2 [21] place significant limits on slepton masses: rngj > 99.0 GeV, 
"^Ai > 91.0 GeV, ruf^ > 85.0 GeV (these assume that the slepton is not nearly-degenerate with 
the LSP) and mi> > 43.7 GeV (from studies at the Z pole). Furthermore, the sneutrino masses 
are closely tied to the left soft mass inputs, and, to avoid extra controversial assumptions, we 
will restrict ourselves to regions of the MSSM parameter space where the LSP is the lightest 
neutralino rather than a sneutrino. To optimise the ino leptonic BRs without running afoul of 
the LEP2 limits, it is best to set m^^ = m^^ . If all three generations have the same soft inputs 
(with Ar = 0), then the slepton sector is effectively reduced to one optimal input value (which 
we identify with fni^). However, since ino decays to tau-leptons are generally not anywhere 
near as beneficial as are ino decays to electrons or muons, it would be even better if the stau 
inputs were significantly above those of the first two generations. This would enhance the 
inos' BRs into electrons and muons. In the general MSSM, we are of course free to choose 
the inputs as such. Doing so would also weaken restrictions from LEP2, especially for high 
tan P values. If we set the soft stau mass inputs 100 GeV above those of the other sleptons 
(with still kept at zero), the lowest allowable slepton masses, presented in the M2 vs. /x 
plane for tan/? = 10 and 20, are as shown in the upper pair of plots in Fig. 1, while if all three 
generations have the same soft inputs we obtain the lower pair of plots in Fig. 1. 

Incorporating such optimal slepton inputs and then scanning over the ino parameters M2 
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Figure 2: BR(iJ^ — > 3iN), where £ = or /i^ and N represents invisible final state 
particles, with (m^, tan/5) = (a) (300 GeV, 10), (b) (500 GeV, 10), (c) (300 GeV, 20), (d) 
(500 GeV, 20). Slepton mass inputs are optimised as in the upper plots of Fig. 1. The shaded 
areas are excluded by LEP; rrit = 175 GeV, rub = 4.25 GeV. 

and /i, for a couple of values of tan/? and m^, yields Fig. 2 for BR(if^ 3iN), where i 
may be either or fi^ and N represents any number of undetectable final state particles 
(either LSPs and/or neutrinos). In these plots, and in plots to be shown hereafter, all possible 
charged Higgs boson decay modes which can result in a final state with three charged leptons 
and no hadronic activity are included, except for leptons coming from tau decays. In this 
figure, including is from decaying taus would not noticeably affect the BRs since the staus 
which could greatly enhance tau production are pushed up in mass. 

As expected, BRs are larger for the rriA = 500 GeV plots on the right than for the = 
300 GeV plots on the left since more 3£-producing H"^ — > inos decay modes open up as 
increases. BRs also decline as tan/5 is raised from 10 to 20. If instead the three slepton 
generations have degenerate soft mass inputs, then Fig. 3 is obtained. Here tua is fixed at 
500 GeV and the left- and right-hand plots depict, respectively, BRs without and with the 
inclusion of is from tau decays. Overall rates drop relative to those in Fig. 2 since: (i) the 
slepton mass inputs must be set higher to evade LEP2 constraints; and (ii) charged Higgs 
boson decays leading to staus via inos - which are now very significant - often result in 
hadronic final states rather than purely leptonic ones. 

Note from Figs. 2 & 3 that low values for are strongly favoured. This can be understood 
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Figure 3: BR(if^ — MN), where £ = or /x^ and N represents invisible final state particles; 
^s resulting from tau decays (are not) are included for plots on the (left-) right-hand side; 
rriA = 500 GeV and tan/5 = 10 in plots (a) and (b), 20 in plots (c) and (d). Slepton mass 
inputs for all three generations are optimised as in the lower plots of Fig. 1. The shaded areas 
are excluded by LEP; rrit = 175 GeV, rrib = 4.25 GeV. 

by inspecting the tree-level decay width formula for — > xfx^ [10], 



T(H 



± 



Fl = cosP[NjiVa + ^,{Nj2 + Njitan9w)Vi2], 

Fr = sinp[Nj^Uii- ^,{Nj2 + Njiiajiew)Ui2l (2) 

where is the sign for the neutralino mass eigenstate (needed to generate a positive physical 
mass given the form of the neutralino mixing matrix), g is the SU{2)j^ coupling and A = 
(m^± — m~± — m~o)^ — 4m~±m~o. Vn & Un {Vi2 & Ui2) give the gaugino (Higgsino) component 

of chargino xf while A^^i & Nj2 {Nj^ & Nj^) give the gaugino (Higgsino) components of 
neutralino x^. We immediately see from F^ and Fr that if the chargino and the neutralino 
are both pure gauginos (the SUSY counterparts of charged Higgs bosons decay into two gauge 
bosons — for which there is no coupling at tree level) or both pure Higgsinos, then the tree- 
level decay width is zero. Simple phase space considerations favouring decays to lighter inos 
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then disfavour situations in which ^ M2 (or <^ M2) in which case hght charginos and 
hght neutrahnos are almost pure gauginos (Higgsinos) — 1/^1^ M2 is preferred, ideally with 
both values as small as possible to make the hghter inos as hght as possible (to the extent that 
LEP2 constraints permit). Thus the optimal region for high inos BRs is where inos 

are mixtures of gauginos and Higgsinos just above the bends of the LEP2 parameter space 
bounds (shaded regions in Figs. 1-3) in the M2 vs. ji plane''. 

In addition, one would like to optimise — > xfx^j decays where j 7^ 1 to obtain the vast 
majority of the decays generating three leptons. Since Mi ~ |-^2, M2 ^ generates an LSP 
that is mostly a U(1)y bino and a somewhat gaugino- dominated chargino — which is bad for 
BR{H^ ^ Xf X?) 

— but also makes for a quite light LSP, which over- compensates for the 
sub-optimal coupling. To increase the other if^ light inos BRs, the mass of the LSP may 
be raised by making M2 somewhat larger than Thus the final pcrscription for optimal 
rates is for small |//| values and shghtly larger, but still small to moderate values for M2. 

The charged Higgs boson BRs must now be tied to the production rate to obtain an 
expected number of signal events. Lowest order (LO) results from the parton-level process 
gb — > tH~ are strongly dependent on which 6-quark Parton Distribution Function (PDF) is 
chosen for convolution and on the scale at which ag is evaluated. Moreover, the 6-quark in 
the initial state originates with a gluon splitting into a bb pair inside the proton, so that 
the above 2 — s> 2 process (when convoluted with initial state radiation involving g ^ bb in 
the backward evolution) can alternatively be taken as the 2 — > 3 hard scattering subprocess 
gg — > btH~ interfaced to gluon PDFs. The two descriptions have complementary strengths: 
the former most aptly describes 'inclusive' tH~X final states, as it re-sums to all orders large 
terms of the form cts log(Q/m;,) (typically Q mt + Mh±, the choice we adopt here), which 
arc absorbed in the phenomenological PDF of the initial 6-quark, while the latter modelling is 
better at describing 'exclusive' observables, as it accounts for the correct kinematic behaviour 
at large transverse momentum of the additional (or spectator) 6-quark in the final state. Yet 
contributions from the two processes cannot simply be summed. In fact, the first term of the 
6-quark PDF is given by the perturbative solution to the DGLAP equation 

6'(x, Q) = ^ log (^) /' ^P,, f ^) g{y, Q), (3) 

where Pgb{z) = {z^ + (1 — ^)^)/2 is the gluon-to-6 splitting function, and the resulting contri- 
bution to gb — > tH~ is already accounted for by gg — > btH~ in the collinear limit. Thus, when 
combining the 2 — > 2 and 2^3 processes, the above contribution should be subtracted from 
the former to avoid double counting [22]. An alternative approach [23] involves specifying 
a threshold in the transverse momentum of the spectator 6-quark, p^""*'^^; then utilising 
2 — > 2 kinematics when p^^ < p'^^^"^ and 2 — > 3 kinematics for p^ > p'^^^^. This is particu- 
larly well-suited to Monte Carlo (MC) event simulations since it does not involve making the 

''For higher values of tan/3, Fl oc cos/3 is small compared to Fr oc sin/?. So the to 5J7(2)^-wino 
Higgsino decay SUSY-related to W^h (where h is now mostly from the down-coupling Higgs doublet 

and so the corresponding Higgsino has a dominating Nj4 component entering into Fl) is also small. But the 
actual inos may not have such compositions. Furthermore, the signature of — > inos is more distinctive 
than that of hW^ — even if the BRs for the two processes were similar, more events from the former 

than from the latter would remain after sufficient cuts were made to eliminate backgrounds. 
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aforementioned subtraction with its associated negative weights. Both techniques yield cross 
section values midway between the larger predictions from gb tH~ and the smaller ones 
from gg — > htH~ (the latter being as much as a factor of 3-4 below the former) . 

Both approaches are less sensitive to the choice of the fe-quark PDF and the mass factori- 
sation scale, Q, than if the two processes were considered separately. However, in each case, 
only some parts of the NLO corrections are accounted for, finally yielding a negative NLO 
contribution. Quite importantly, recent results [24] have proved that full NLO corrections to 
the 2 — > 2 process (i.e., including both one-loop and radiative QCD corrections) can yield 
an overall /C-factor much larger than one^, overturning the negative corrections obtained via 
the above procedures. Thus it is no longer justifiable to adopt normalisations based on these 
techniques. 

The MSSM implementation [25] of the HERWIG [26] event generator was used to simulate 
the gh — > tH~ process and the various backgrounds. As most backgrounds are only known to 
LO accuracy, no additional /^-factors were incorporated and default LO PDFs and Ug values 
were employed. This partly explains the improvement to be seen herein relative to Ref. [3], 
where normalisation was via the old subtraction procedure. Nonetheless, we still regard our 
results as conservative since the dominant backgrounds (after cuts) are tt production, which 
has a similar QCD X-factor to that of the signal, and irreducible contributions from direct 
neutralino-chargino pair production, which, being EW processes at tree level, have smaller 
QCD corrections (of the order of 20% or so [27]). Yet one should also verify that the additional 
6-quark at high transverse momentum produced by the gg — > htH~ contribution, which is not 
present in the (infrared dominated) backward evolution of the 2 — > 2 process' initial fo-quark, 
does not render untrustworthy a kinematical analysis done solely utilising the 2 — > 2 process. 
We have confirmed this by also running HERWIC with gg — ^ htH~ as the hard subprocess, 
adopting our usual selection cuts, and checking that in fact observable quantities (distributions 
and event rates) are not significantly afi'ected by the presence of a spectator 6-quark in the 
detector. All results shown will correspond to the outputs of the 2 — > 2 process. 

Fig. 4 shows a{pp tH~X + c.c.) with the subsequent decay — > 3iN, where M2 and 
H are fixed at the favourable values of 210 CeV and 135 GeV, respectively — leading to the 
exclusion of tan/5 values below ~ 5 by the 103 CeV LEP2 lower bound on the chargino mass 
[21]^. In the plot, the preference for high and low values of tan/3 so well-known for the raw 
pp tH~X cross section remains, though rates are nevertheless sufficient to seek a visible 
signal even in the intermediate tan (3 region via our characteristic signature. 

It is instructive to next isolate the dependence of the signal rate upon the masses of the 
sleptons^°. This is done in Fig. 5 for three choices of the other relevant MSSM parameters. 

^The size of the X-factor is sensitive to the choice of the mass factorisation scale, Q, and the couphng 
renormaUsation scale. Equating these two scales for convenience, -ftT-factors near unity are found for low scale 
choices, Q ^ ^{mt + mH±), while higher scale choices can yield values on the order of 1.6-1.8. See [24] for 
further details. In all cases a K-iactoi much less than one is not found. 

^For this choice of input parameters, the m~± bound is probably more restrictive than the one from 

Higgstrahlung, e+e~ hZ (and hA); however, this will not be true for other choices of M2 and /i, such as 
those considered in the next paragraph. Note that the location of the Higgstrahlung bound is quite vague due 
to uncertainties in the radiatively-corrected mass mh and errors in the measured value of mj. 

^°Though sleptons are light, direct BRs to slepton pairs are at the sub-percent level. Sleptons mean- 
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Figure 4: a{j)p tR-X + c.c.) x BR(i7± ^ UN) (in fb), where I = or and 
represents invisible final state particles: M2 and n arc as noted, and £s from tau decays are 
included. Slepton mass inputs are optimised as in the upper plots of Fig. 1. The LEP2 M~± 

Al 

hmit excludes the shaded region, and the on the left signifies where the BR is virtually zero. 

All combinations fix tan/? at 10 and at 500 GcV. Parameter Set A (solid curve in Fig. 5) 
also sets n = 135 GeV and M2 = 210 GeV (as in Fig. 4), while Parameter Set B (thick dashed 
curve in Fig. 5) has = 150 GeV and M2 = 280 GeV and Parameter Set C (dot-dashed curve 
in Fig. 5) adopts /x = —150 GeV and M2 = 300 GcV (these same parameter sets will also 
be used in the forthcoming experimental analysis). The horizontal axis in Fig. 5 is the soft 
slepton mass input (as before, left and right soft masses are degenerate and A-terms are zero). 
Bear in mind that this is not the same as the physical masses of the various sleptons, which 
also have so-called D-term contributions. The curves are terminated on the left side at the 
point where LEP experiments preclude the resulting light sleptons. Also shown by the dotted 
curve is the effect of removing the equality m^^ = fni^: in this case m^^ = m^^ -|- 100 GeV 
while all the other MSSM parameters arc the same as in Parameter Set A. 

Focusing the account upon the two curves relating to Parameter Set A (features of the other 
curves are seen to be qualitatively similar), a sharp drop is seen around m^^ ~ 123 — 125 GeV 
where the second neutralino becomes degenerate with the charged selectrons and smuons 
and also where the lighter chargino becomes degenerate with the sneutrinos (of the first two 
generations). The drop is due to the closing of the two-body decay modes X2 i^^i^ and 
Xi ~^ 1 where £^ and vg, are on-mass-shell. Although the two-body decay modes close at 
this point, the sleptons still make their presence felt in the associated three-body decay modes 

ingfuUy influence charged Higgs boson leptonic BRs via the sleptons' involvement in subsequent ino decays. 
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Figure 5: a{j)p tR-X + c.c.) x BR(//± ^ UN) (in fb), where £ = e± or ji and repre- 
sents invisible final state particles, vs. rrig^, the soft slepton mass input for the first two genera- 
tions (soft stau mass inputs are pushed up by an additional 100 GeV, Ar = 0). The set (M2, /i) 
is fixed at (210 GeV, +135 GeV) for Parameter Set A (solid curve), at (280 GeV, +150 GeV) 
for Parameter Set B (thick dashed curve) and at (300 GeV, —150 GeV) for Parameter Set 
C (dot-dashed curve). The dotted curve replaces m^^ = m^^ in Parameter Set A with 
m^^ = + 100 GcV. The curves are terminated at the left where they would be LEP2 
excluded, including the additional condition that mc > m~o; is from tau decays are included. 

via off-shell contributions. A modest rise in the rate occurs for the solid curve as reaches 



~146 GeV where the second neutralino and the sneutrinos (of the first two generations) become 
degenerate and the 'spoiler' modes, X2 ~^ ^i^i — which result in no charged leptons — become 
inaccessible, consequently allowing BRs for the ^-producing channels to rise. This feature is 
absent for the dotted curve, with m^^ = m^^ + 100 GeV, since two-body ino decay modes to 
the now too heavy sneutrinos are not open. Yet despite the absence of the spoiler modes, rates 
remain lower in this case because i^s are also heavy, simultaneously weakening the rate to 
charged leptons (when both are accessible and neither is phase-space suppressed, ino decays to 
either £ls or to i>£S may be larger, depending on the composition of the neutralino). Note also 
that the dip at m^^ ~ 123 — 125 GeV is less pronounced since now only two-body decay modes 
to IrS are turning off at this point rather than to both left and right charged sleptons as in the 
solid curve. Lower values of m^^ are now possible since our (perhaps unnecessarily restrictive) 
requirement that > m~o is satisfied for lower m^^ values. However, even going to such 
modest m« values does not compensate for the enhancement obtained when = m« 



and the two-body modes to left charged sleptons arc available. Thus the peak magnitude is 
appreciably lower for the dotted curve (confirming that = mi is the optimal setup) 
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Setting > m^^ shifts the curve to the right and sUghtly lowers the peak plateau. 
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For Parameter Set A, with set to 110 GeV, the largest contributor to the signal 
events is in fact^^ xtx2 (35!2%), followed closely by xtxl (34.8%) and then 

xtxl (19.3%), xtxl (6.5%), and small contributions from xfxl (3-0%) 

and iJ^ xtxi (1-2%) (here, the three is all come from the xf)', — > xfx2 3 channels do 
not lead to most of the prospective signal events, contrary to what was assumed in Ref. [2]. 
Nevertheless, rates seen in Fig. 5 arc still closely linked to m~± and m~o since Xa and/or xH 
are present in most (92.3% for Set A) events, and, even if one (or both) is not in the ino pair 
to which directly decays, the heavier inos into which does decay in turn sometimes 
decay into these lighter inos (and charged leptons or neutrinos) to generate the signal events. 
With so many contributing channels, some of which involve multiple sparticle to sparticle 
decay chains, simulation of the signal with a robust event generator is imperative to ascertain 
the percentage of the events predicted utilising sparticle BR assignments that survive the cuts 
needed to sufficiently identify the signature and eliminate the backgrounds. 

Returning now to the question of potential backgrounds from colourcd-sparticle production 
processes, gluinos and squarks of the first two generations may in principle produce multi- 
lepton events with top quarks; however, in practice, top quarks are quite often not present 
in such events. Further, the limit on the squark (gluino) masses from Tevatron studies is 
now at least ~260GeV (~ 190 GeV) [28], and will rise if Tevatron searches continue to be 
unsuccessful. In addition, if the gaugino unification assumption also encompasses the gluino, 
then the gluino mass would be in the range ~700-1000 GeV for the points being considered, 
and (at least in mSUGRA-inspired scenarios) squarks are expected to have heavier or at least 
comparable masses [29]. Thus there is substantial rationale for limiting this analysis to heavy 
gluino and squarks (of the first two generations) masses. 

Stops are different though. Stringent experimental limits from LEP2 on stop masses only 
set a lower bound of ~100GeV [21], and stop pair production will generally lead to events 
containing top quarks. Possible decay chains that could mimic our signal events include for 
example U* Xiit ^ biu) + {xl ^ X?^^)(^ ^ hadrons) and U* {xt Xi^b + 
(X2 ~^ Xi^^){'t hadrons). Sbottom pair production can also yield such a final state, as 
for example via bb* — > {xt — > Xi'^^){i ~^ hadrons) + (x2 Xi^^)b- Note though that 
such processes do have an extra 6-jet (typically with high Pt) beyond that expected from 
gb {t ^ hadrons) {H~ — > x7X^ ~^ 3iN) where N may be any number of colourless neutral 
stable particles. Fortunately, our studies indicate that the extra 6-jet that is present in the 
2 — > 3 charged Higgs boson production process tends to be rather soft. So a cut on extra hard 
jets in the event does tend to remove the background from stop and sbottom pair production 
(as well as that from squark and gluino production in general). 

In keeping with the optimal strategy outlined above for the slepton sector, stops are made 
heavy to minimise this potential background. Thus we deal only with the MSSM backgrounds 
that must be present: that from direct ino pair production (since we require — > inos, the 
inos must be relatively light), and what coloured-sparticle backgrounds still remain when we 
have heavy gluinos and all squark inputs pushed up to 1 TeV. A more in depth study of light 
stops possibly mimicking our signal will be presented in an upcoming analysis [30]. 

^^Numbers include leptons from decaying taus, but said inclusion only causes slight changes. 
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3 mSUGRA Parameter Space 



Before initiating the experimental analysis, we would like to document the potential for util- 
ising the '3£ + t ' signature from H"^ inos in the more restrictive mSUGRA parameter 
space. As we will soon see, here prospects are quite bleak. In mSUGRA, the free param- 
eters are generally set as tan /3, a universal gaugino mass defined at the Grand Unification 
Theory (GUT) scale (-^i), a universal GUT- level scalar mass (Mq), a universal GUT-level 
trilinear scalar mass term (Aq), and the sign of As already noted, the signal has a strong 
preference for low values of Yet in the mSUGRA scenario, is not a free parameter 
— it is closely tied to the masses of the scalar Higgs bosons via the Mq input. Furthermore, 
the different soft slepton mass inputs can no longer be set independently: in particular, when 
evolved down to the EW scale using renormalisation group equations, the staus' soft inputs 
tend to be lower than those of sleptons from the first two generations rather than higher as 
was put in by hand in the more favourable MSSM parameter set choices of the preceding 
section. Fig. 6 shows the values for a{pp tH~X + c.c.) x BR(i/^ — > "ilN) obtained for 
several discrete values of tan /3 and > (analogous plots for /i, < are similar) with Aq 
set to zero. The excluded regions shown take into account constraints from LEP2 save that 
coming from Higgstrahlung^^, but not additional constraints^^ from b ^ s^, — 2 and other 
loop-level effects (nor considerations from cosmology) which are now harder to dismiss since 
the behaviour of the model is specified all the way up to the GUT scale. 

In regions of parameter space not LEP2-excluded, maximal rates are found for very high 
values of tan (3. Here, charged Higgs boson decays into three leptons are most favourable when 
Mo is low and Mi is high. However, the production process (which grows with tanP once 
the minimum around tan P — 6 is passed) favours low Mq and low Mi . Therefore there is an 
optimal Mq value for the maximal signal rate (in regions not excluded by the afore-mentioned 
LEP2 constraints) which grows from ~300 GeV to ~500 GeV as tan f3 is increased from 20 to 
40. Even larger rates for yet higher values of tan/? (>60) are swallowed up by the greatly- 
expanded LEP2-excluded region (sparticles or Higgs bosons become unacceptably light and/or 
the fighter stau becomes the LSP) . Inputs must be chosen to make soft stau masses high (while 
one would like to — but in mSUGRA cannot — keep the other soft slepton inputs low to 
obtain good leptonic decay rates) to avoid a stau LSP and the LEP2 bound on the physical 
stau mass. Thus the window of allowed points for such extremely high values of tan/3 is 
orthogonal to where substantial signal rates are possible. 

Even with 500 fb~^ of integrated luminosity, the handful of signal events expected in the 
best allowed cases would probably be unresolvable from amongst the backgrounds. Thus, 
given the very meager chances of extracting a signal even in the perhaps overly-generous 
allowed regions of parameter shown here, a more thorough mSUGRA analysis would probably 
be irrelevant. 

^^It should be noted that the small unexcluded regions shown where the cross section is ~ 0.01 fb may 
be partly or totally excluded by the LEP2 Higgstrahlung constraint. However, given the aforementioned 
uncertainties surrounding this limit, we conservatively make no attempt to place corresponding exclusion 
contours on our plots. 

^^See [31] for a more complete analysis of the present-day constraints on the mSUGRA parameter space. 
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Figure 6: a{pp tH^X + c.c.) (in fb) multiplied by BR(i/^ 3iN), where i = or fi^ 
and N represents invisible final state particles, for a spread of mSUGRA parameter sets in 
the Mq vs. Ml plane; = in all plots and £s resulting from tau decays arc included. The 
number in the upper right of each plot is the tan f3 value followed by the sign of fi. The solid, 
dot-dashed, and dotted contours are for 10^^ fb (also labeled). 10^'^ fh and 10^'^ fh, respectively. 
The shaded regions are excluded by theoretical considerations or LEP2 measurements (save 
that constraints from Higgstrahlung are not applied). 



4 Experimental Analysis 

In the previous sections we outlined the potential for observing the charged Higgs bosons 
through their decays into charginos and neutralinos, eventually yielding three leptons plus 
missing energy, and in the presence of a hadronically reconstructed top (anti)quark. As a 
next step, we study the feasibility of detecting such a signal in a realistic LHC detector 
environment (CMS). We use the MC event generator HERWIG (version 6.3) and simulate the 
gb — * tH~ + c.c. — > 3£ + p™'^^-|-t signal for the three MSSM settings already discussed, which 
we specify more fully here: 



Set A: M2 = 210 GeV, n 



135 GeV, m^^ 



110 GeV, m~g = 800 GeV, m„- = 1 TeV. 
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• Set B: M2 = 280 GeV, = 150 GeV, m^^ = 130 GeV, Trig = 900 GeV, rriq = 1 TeV. 

• Set C: M2 = 300 GeV, // = -150 GeV, m^^ = 150 GeV, m-g = 1 TeV, rrig = 1 TeV. 

Recall that in all settings we assume Mi = |tan^^vK-^2- Furthermore, for sleptons and 
squarks we will always take soft mass inputs for all generations to be degenerate (with m^^ = 
m^^). The physical sneutrino masses, m^, can be derived from the above parameters and are 
approximately 90, 115 and 135 GeV for the respective scenarios (when tan/3>5). Parameter 
Set A lies inside the optimal region in the three-dimensional (M2, /i, ttii^) space identified in 
Sect. 2, whereas Set B is a more borderhne case and Set C is a difficult case with a negative ji 
parameter. Set A features light inos and sleptons, allowing several supersymmetric decay 
modes to have considerable BRs for relatively moderate values of itla (and mH±)- The ino 
sectors in Set B and Set C are heavier, thereby limiting the number of possible sparticle decay 
modes. In Set B sleptons are light, whereas in Set C these sparticles are also heavy. This 
last difference markedly alters the kinematics in ways to be discussed shortly. The MSSM 
sparticle spectrum and decays are obtained from ISASUSY 7.58 [32] through the ISAWIG 
interface [33]. ISASUSY contains a one-loop treatment of all Higgs boson masses and tree- 
level sfermion masses. Several three-body decays are included, taking into account the full 
Yukawa contributions, which are important in the large tan (3 regime. The charged Higgs 
boson BRs are taken from HDECAY [34] (again, via the ISAWIG interface), which calculates 
these in accordance with the most recent theoretical knowledge. For the SM backgrounds, all 
leading processes that can produce the 3£-|- p™*^*^ + t signature have been simulated: tt {thW~ 
typically is 1/4 as large), ttZ, tt7* and flh. Furthermore, all SUSY backgrounds have been 
considered for the chosen settings: ino pair production (including squark-|-ino production), 
squark and/or gluino production and slepton pair production. Of these, the first listed class 
of SUSY contributions has the largest cross sections in general, because inos are fairly light in 
comparison to the coloured sparticles. In our scenarios, slepton pair production never results in 
a three lepton final state^"'', so that it will be excluded from further consideration. The detector 
aspects were simulated using CMSJET 4.801 [35], which contains fast paranietrisations of the 
CMS detector response and, for 6-tagging, a parametrised track reconstruction performance 
based on GEANT simulations [36]. 

In Parameter Set A, the neutralinos and have masses of 78, 131, 146 and 

253 GeV, respectively. The masses of the charginos Xi ^iid xt ^'^^ and 252 GeV. The 
are allowed to decay into all kinematically accessible ino pairs, xfx°, which in turn 
can decay into three leptons^^ (electrons and/or muons) plus invisible neutral particles (xiS 
and/or neutrinos). In this scenario, the primary source of three-lepton events (before any 
kinematical cuts are considered) is generally charged Higgs boson decays to xfyH- This is 
true for 225 GeV < 400 GeV and^'' tan/? ~ 10. Charged Higgs boson decays to xt^H 
and Xi^A ^-^s also important sources of 3£ events in this region of parameter space, and 
the contributions from these modes grow to equal or surpass that from X2 decays for 
ruA ~ 400 GeV. X2 decays almost exclusively via an intermediate state containing an on-shell 

^^Unless four leptons are produced rather than the usual two, and then one lepton is subsequently disregarded 

due to having a pT vahie too low to pass our cuts. Rates for such events are negligibly small. 

^®Here i = 1,2 j = 1,2,3,4; if j — 1, then the three leptons must all come from cascade decays of the 
chargino. 

^''The upper (lower) value drops by ~20GeV (-^lOGeV) as tan/3 goes from 10 to 30 (5). 
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charged slepton, while xt decays through a intermediate state including an on-shell sneutrino. 
Here though BRs for stau decays rise as tan P grows (in part due to the fact that, for fixed 
soft slepton mass inputs, the physical mass of the decreases swiftly as tan /3 rises to higher 
values): for tan (3 = 5, 10, 20, 30, the X2 BR to staus is about 0.35, 0.42, 0.61, 0.74, respectively. 
Additionally, about 1/3 of the xt decays to sneutrinos are to z/^, and Xi decays to rf^z^r become 
accessible at tan/? ~ 11 with the BR for this decay growing to 0.24 (0.06) when tan/? reaches 
30 (20). These ino to stau and V^- decays reduce the number of X2 ~^ -^^-^^ and Xi ~^ ^d"^ 
decays, where £ = e or //, which lead to virtually all the signal events that survive the necessary 
cuts"*^^. Leptonic tau decays arc allowed in the event generation, although daughter leptons 
from these will mostly be rejected during the analysis stage due to their softness (low pr^). 
Crucial mass differences have values^^ of (m~o —mj^, iti~± —m~ , mj^ —m~o) = (~10-15 GcV, 

'*-'l ^1 

~10-22GeV, ~35-45GeV). In all these cases there is enough phase space for most of the 
resulting leptons to have sufficiently high transverse momenta. 

In order to distinguish between the signal and the backgrounds (both SM and MSSM), 
we will apply a set of selection criteria that will allow us to obtain a favourable signal-to- 
background ratio using only physically well-motivated cuts {i.e., with only a very loose de- 
pendence upon the MSSM parameters). We will first explain the selection strategy and then 
illustrate the results numerically in a table. 

First of all we require the following basic topology: 

• Events must have exactly three isolated leptons {i = e,//) with pt > 20, 7, 7GeV, all 
with |r/| < 2.4. The isolation cut demands that there are no charged particles with 
Pt > 1.5 GcV in a cone of radius AR = ^ (A^)^ + (At^)^ = 0.3 radians around each 
lepton track and that the sum of the transverse energy deposited in the electromagnetic 
calorimeter between AR = 0.05 and AR = 0.3 radians be smaller than 3 GcV. 

The choice of the minimum p^ value for the leptons is driven by both trigger and back- 
ground rejection considerations. In the case of muons, requiring a hardest lepton above 20 GeV 
is already sufficient for the event to be triggered with 90% efficiency by the single-muon trigger 
under low luminosity running conditions at the LHC (for electrons this threshold is somewhat 
higher) [37]. Apart from the single leptons triggers, the di- and trilepton thresholds will in- 
crease the efficiency for triggering on the 3i signal. The tight isolation criterion is needed in 
order to reject leptons coming from heavy flavour decays, especially in the low pr region. As 
we will discuss later, it is very effective against, for instance, the tt background, when one or 
more of the leptons originates from a 6-jet. 

Apart from requiring the three leptons, it is also necessary to reconstruct the (hadronically 
decaying) top quark that is produced in association with the boson. This is mainly 
motivated by the need to strongly suppress the tt and ino-ino backgrounds. A reconstructed 
top (antitop) quark is recognised via the following cuts: 

• Events must have at least three jets, each with px > 20 GeV in I77I < 4.5. 

^^If the soft stau mass inputs are made heavier, rather than degenerate with the other soft slepton mass 
inputs as is done in this analysis, the unprofitable stau channels could be eliminated and the number of events 
could as much as double for high values of tan /3. 

^^These numbers depend moderately on tan p. Values given here (and later for Sets B & C) cover the range 
of interest in this work: 5 < tan /? < 30. Physical slepton masses are those of the first two generations. 
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• Among these, the three jets that are most hkely to come from a top quark decay are 
selected by minimising rrijjj — rrit, where rrijjj is the invariant mass of the three-jet 
system. This invariant mass rrijjj must be in the range rrit ± 35 GeV. 

• Two of these three jets are then further selected by minimising mjj — Mw±. Their 
invariant mass, rrijj, must be in the range Mw± ± 15 GeV. 

• The third jet {i.e., aside from the two jets in the preceding point) must be 6-tagged. 
(This we equate with a significance of the transverse impact parameter a{ip) = 
which is larger than 2.) 

A strong rejection of tt events is obtained after the requirement of a hadronically recon- 
structed top quark in addition to the three leptons. Assuming that the three jets recon- 
structing rrit a-re indeed correctly assigned, this requirement means that the second top should 
provide two leptons (one from the and one from the b) while the third lepton should come 
from initial/final state radiation {b,K,7i, ...). In this case, two leptons will be in general soft 
(< 5 GeV) and non-isolated. Another scenario in which ti production can lead to a 3i + t 
final state is the one where both top quarks have decayed leptonically, and two radiated jets 
accidentally reconstruct the mass and then combine with a 6-jet from top decay to mimic 
a hadronically-decaying top quark. Here, two leptons can be hard, but the third one must 
still be soft and in general non-isolated. Therefore, in order to achieve a sufficient suppression 
of the ti background, we have chosen to set the lower limit on the pt of the leptons at 7 GeV 
(although lowering it would increase the signal yield) and to apply a tight isolation criterion. 

Whereas the ti background is greatly suppressed by the previous selection steps, tiZ, tiY 
and tih events would still survive the 3i+t criteria. Therefore we require an additional Z-\eto: 

• Reject all events with di-lepton pairs with opposite charges and the same flavour that 
have an invariant mass in the range Mz ±10 GeV. 

The Z-veto rejects tiZ events efficiently. Moreover, although the ^7* and tih backgrounds 
largely survive this requirement, their residual cross sections are now innocuously small. 

In addition to eliminating the SM noise, cuts to efficiently suppress the SUSY backgrounds 
that can lead to a 3i + t final state must be considered. As mentioned before, slepton pair 
production does not pose a problem in our scenarios since it cannot lead to a three- lepton final 
state. Ino pair production and squark-|-ino production can have large cross sections; however, 
most events from these processes do not contain a top quark and will thus be rejected via the 
hadronic top requirement. Events that are still left after this cut form the main irreducible 
SUSY background. Squark/gluino production is another potentially dangerous source of noise. 
These events, however, typically contain many energetic jets besides those coming from the 
top decay (as previously intimated). Therefore, they can be rejected using an additional jet 
veto: 

• Reject all events containing any jets (other than the three jets selected for the top 
reconstruction) with pt > 70 GeV and \r]\ < 4.5. 

For further signal-to-background rejection, we impose the following (slightly model depen- 
dent) selection criteria (here optimised for rriA — 350 GeV and tan/3 = 10): 
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• For the three isolated leptons aheady selected, the pt of the hardest lepton should be 
below 150 GeV whereas the pr of the softest lepton should be below 40 GeV. 

• The missing transverse energy should be larger than 40 GeV. 

• The effective mass, M^g, constructed from the p^ and p™'^^ vectors as M^g = 
^J2p^p™^^^{l — cos A0), is required to be lower than 150 GeV (here is the azimuthal 
angle between pj: and p™'®®) . 

The missing energy requirement has little affect on the signal yield since the two x\ 
the final state usually supply sufficient p^^^^] on the other hand, this cut does reduce the SM 
ttV {V=Z,Y) backgrounds. As was shown in [2], the effective mass variable does have some 
dependence on the ino mass spectrum; but it also proves to be effective against the above ttV 
processes plus squark/gluino and ino pair production backgrounds as well. 

After applying these selection criteria, we obtain the number of signal (S) and background 
{B) events given in Tab. 1, assuming Parameter Set A, with rriA = 350 GeV and tan/3 = 10, 
and for an integrated luminosity of 100 fb~^. Results shown therein clearly confirm the points 
made in the preceding description of the cuts. 



Process 


3£ events 


Z-veto 


hadr. top^ 


6-tag 


jet veto^'^ 


others''"''^ 


ti 


2781 


2465 


91 


15.5 


11.1 


5.8 


ttz 


492 


82 


19 


8 


2.4 


0.8 


tij* 


22 


21 


7 


2 


0.4 


0.2 


tth 


59 


52 


17 


4 


1.6 


0.2 


XX 


19993 


18880 


237 


31 


9 


3 


g 


12712 


11269 


3984 


861 


6 


1 


tH- & tH+ 


508 


485 


126 


36 


29 


25 



Table 1: Number of signal and background events assuming Parameter Set A, with = 
350 GeV and tan/3 = 10, for 100 fb~^. (Note that the difference between event rates in the 
'hadr. top' and the '6-tag' columns is not only due to the experimental 6-tagging efficiency 
but also takes into account part of the algorithmic efficiency.) '''Here, top reconstruction 
requires > 3 jets with pt > 20 GeV, rrijj ~ My^^± and rrijjj ~ . ''"'"Here, one vetoes 
additional jets beyond 3 with px > 70 GeV. '^''"^Here, one imposes Pt(^i) < 150 GeV, prih) 
< 40 GeV, pf^^ > 40 GeV and M^s < 150 GeV. 

Fig. 7 shows the three-lepton invariant mass distribution for our typical signal (m^ = 
350 GeV and tan/3 — 10) on top of the background (SM + SUSY) for an integrated lumi- 
nosity of 100 fb~^. The peak in the three-lepton invariant mass distribution depends both on 
mH± Bind on the mass spectrum of the intermediate charginos and ncutralinos. Therefore, 
a direct 'parameter-independent' mass reconstruction does not seem feasible at this stage. 
The determination of the charged Higgs mass will require comparisons between the measured 
three-lepton invariant mass and MC distributions. 

Especially since we cannot claim ability to discern a mass resonance, it is well to review the 
strategy of this work. We attempt to locate a signal for a charged Higgs boson by comparing 
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Figure 7: Three-lepton invariant mass distribution for Parameter Set A, with niA = 350 GeV 
and tan/3 = 10. The signal peak is shown on top of the SM + SUSY background for an 
integrated luminosity of 100 fb~^, after all cuts described in the text. 



the number of events fitting our criteria at a specific point in the MSSM parameter space with 
and without inclusion of our gh tH^ (and c.c.) subprocess accompanied by if^ — > ?>iN . 
Events resulting from sparticle production processes (sleptons, squarks, gluinos, and/or inos) 
which happen to satisfy our criteria are regarded as backgrounds. It is conceivable that the 
rate from sparticle production processes at a different point in the MSSM parameter space 
could mimic the excess due to charged Higgs boson production at the afore-mentioned specific 
point we are considering. Thus one could consider instead including the sparticle production 
processes in the signal category and mapping out the excess expected from all non-SM sources 
throughout the MSSM parameter space. 

We have chosen to proceed treating only '''-events as the signal events for several reasons: 
(1) Our choice of cuts to select signal events is tailored to pick out gh tH~ (and c.c), 
H"^ 3iN events; (2) We do give only cursory consideration (for the moment) to some rather 
limited regions of the parameter space, such as those containing light sbottoms or stops, which 
might yield significant numbers of non-H^ signal events (though we do give reasons why we 
suspect this will not be the case); and, most importantly, (3) we expect that results from 
our channel will be correlated with results in other channels to resolve most of any ambiguity 
about the correct location in the MSSM parameter space. Sparticle production has other 
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signatures. For instance, we cannot restrict our consideration to parameter set choices where 
rates for direct ino production are low since we require the charged Higgs bosons to decay 
into inos. But there are other signatures for direct ino production [38], notably trilepton or 
like-sign dilepton signals (without an associated t or t) from chargino-neutralino production. 
In fact, knowledge from these other channels (either from the LHC or perhaps from runs at 
the Tevatron) could enable future analyses to sharpen the cuts employed here (which only 
incorporate some vague assumed form for the ino spectrum) , and then perhaps reconstruction 
of the charged Higgs boson mass will be possible. Therefore, our perspective is to treat this 
work as one piece of a body of analyses with which we hope to pin down SUSY's nature. 

Maintaining the MSSM setup of Parameter Set A, we can now perform a scan over ttla and 
tan (3 in order to determine the discovery potential for the '3£ + p^^^^ + t' signature we have 
been considering. We assume an integrated luminosity of 100 fb~^ and require the significance 
of the signal, S/Vb, to be larger than 5. The resulting 5(j-discovery potential is shown in 
the top plot of Fig. 8. The left edge of the potential discovery region at ttia ~ 250 GeV is 
determined by the kinematic requirement that mH± > itl-o + m~±. The upper edge in tan/? 
originates from decreasing xfx^ ^ind X? (^^ = Ij 2 and j = 2, 3, 4) BRs. This 

is in part a consequence of the high tan f3 enhancement of couplings to the third generation 
(taus, top and bottom quarks), which grow at the expense of the couplings to the inos (the 
intermediates we need to obtain our hard, isolated electrons and muons) and in part due to 
the increased BRs for ino decays into staus, which grow at the expense of decays into lepton(e 
and /i)-yiclding selectrons and smuons. The upper edge in rriA and lower edge in tan/? are 
determined by the tjia and tan/3 dependence of the production cross section. Conservative 
LEP exclusion limits [39], mainly from Higgstrahlung {i.e., e+e" — > hZ and e'^e" — > hA), are 
also drawn in the figure along with a horizontal dotted line below which m~± does not respect 
the LEP2 bound^o. 

Parameter Set B produces the following mass spectrum: the neutralinos Xi) X2) X3 ^ind 
X4 have masses of 103, 159, 165 and 311 GeV, respectively, while the masses of the charginos 
Xi and xt 131 and 311 GeV. The 5(T-discovery potential for this setting, again for an 
integrated luminosity of 100 fb~^, is presented in the middle plot of Fig. 8, after the usual 
selection procedure. A noticeable difference with respect to Parameter Set A is that here 
the discovery zone starts at somewhat higher values of tua due to the higher tjih^ threshold 
needed for decays to ino pairs, since they are heavier than in the previous scenario. For 
tua ~ 300 GeV, — > xfxs is the dominant source of 3£ events, rather than — > X1X2 
in the previous scenario^^. There are two reasons for this change: firstly, with Set B, xl has 
a more favourable gaugino/Higgsino mixing than does X2'i secondly, for Set B BRs for 
X2 decays into sneutrino spoiler modes for tan/3 = 5, 10, 20, 30 arc about 85%, 84%, 68%, 39%, 
whereas for Set A these values are all roughly 0.5%. H"^ X1X3 decays remain the dominant 

^°Not shown on the plots in Fig. 8 arc upper tan /? bounds of 32.2, 28.9, and 44.0 for Parameter Sets A, B 
and C, respectively, above which the lighter stau mass dips below the LSP (xi) mass. This bound may be 
evaded by raising the soft stau mass inputs above those of the first two generations of sleptons. 

^-'^For Set B, there is a thin strip of parameter space around itla ^ 280-290 GeV in which X2 is the 

dominant source of 3£ events. However, the overall 3£ BR drops precipitously in this region as niA decreases, 
and so there is no potential for discovery. 
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source of 3i events even for high charged Higgs boson masses. The H"^ — > X2X2 ^^'^ ~^ 
XiX\ decay modes, which led to the majority of the three lepton events for a ~450GeV 
charged Higgs boson with Set A, in Set B are no longer dominant for higher m^ji values, at 
most providing ~30% of the events (before cuts) for mH± ~ 650 GeV (for which mass value the 
production rate is already too low for any hope of discovery). The %!] decays predominantly 
into a charged slepton; the BR for X3 decays into staus grows from ~32.5% to ~47.5% as 
tan (3 goes from 5 to 30, cutting into the desired decays to selectrons and smuons. Though Xg 
decay modes to sneutrinos are accessible too, the combined BR for such spoiler modes remain 
at the 3-4% level. The Xi decays through a sneutrino intermediate state and the associated 
charged lepton^^. Crucial mass differences have values of (^j^o —mj±, Tn~± —m~^, mj^ —m~o) 
= (~25-30GcV, -12-22 GeV, ~30-40GeV). As in Set A, there is enough kinematical phase 
space for most of the resulting leptons to have sufficiently high transverse momenta to pass 
the signal selection criteria. 

In Parameter Set C the neutralinos Xv X2! X3 ^^^id Xa have masses of 118, 162, 171 and 
324 GcV, respectively. The masses of the charginos Xi ^-nd arc 143 and 324 GeV. Scanning 
over rriA and tan (3, after the customary selection cuts, now leads to the Sa-discovery potential 
seen in the bottom plot of Fig. 8. The reach both in rriA and tan/3 is strongly reduced in 
comparison to the previous scenarios, in part due to the heavier ino mass spectrum which 
gives the expected upwards shift of the left edge in uia- As with Set B, xfxt is the 

dominant source of signal events for myi<600GeV (m^<520GeV) and tan/? ^ 5 (30). For 
350 GeV < rriA < 450 GeV, virtually all (> 90%) signal events come via this channel. For 
higher masses, — > X1X4 and — > X2X3 contributions grow to become comparable. Also, 
as with Set B, the X3 decays predominantly into a charged slepton; again decays into staus — 
BR ~29% (~58%) for tan P = 5 (30) — cut into the desired decays to selectrons and smuons. 
The sneutrino spoiler modes also have a combined BR roughly in the 10-20% range. The xt 
decays through a sneutrino intermediate state: for tan/3 < 20, about 2/3 of the time into 
sneutrinos of the first two generations and about 1/3 of the time into a Ur and the associated 
T-lepton. For higher values of tan/3, the xf ~^ ^^t^r decay mode becomes accesssible and 
reaches a BR of almost 70% by the time tan /3 reaches 30. Now crucial mass differences have 
values of {m~o —mj^, m~± —'m-^, mj^ —m~o) = (~13-17GeV, ~4-10 GcV, ~30-45 GeV). Most 
significantly, there is considerably less phase space available to leptons produced in chargino 
to sneutrino decays^"^; thus, said leptons are typically too soft and usually fail the pt cut. This 
explains the much smaller discovery reach for Set C compared to the one for Set B. 

For II < 0, the same magnitude of leads to heavier inos (in particular, the LSP and 
lighter chargino). Thus, for a fixed \fi\, we expect a smaller signal rate for fi < than for 
H > 0. However, the more rapid rise of the chargino mass as increases with fi < also 
means that we can go to smaller values on this side before we run afoul of the LEP2 
excluded region^^. Thus, one can shift to lower values on the n < side to obtain roughly 

^^Chargino decays to t^Vt become significant for higlier values of tan/3: this BR is ~21% (~6%) for 
tan /3 = 30 (20). Not so useful decays to i?T-r=^ have BRs of 34-37%. 

^■'Thc difference between m~± — m~ for Set C and the values for Sets A & B is more striking when tan/? 

is restricted to be > 10. Then for A and B m~± — m~ ~17-22GeV while for C the value is ^4-7 GeV, with 

this mass difference growing with increasing tan/3 for A and B and shrinking with increasing tan/3 for C. 
^^This is traceable to a term oc 2/iM2 sin 2/3 in the formula for the chargino mass, and hence the asymmetry 
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the same rates as found on the > side (c/., Fig. 4 of [2]). 

Some perspective as to the new regions of MSSM parameter space that might be probed 
via the ^Si + p^^^ + f channel is provided by Fig. 9, which shows the reach of this — > inos 
signature in the case^^ of Parameter Set A together with those of the H~ r~z/^ and H~ bt 
channels, with tan (3 plotted on a logarithmic scale to better illustrate the intermediate tan f3 
regime. The discovery reaches for channels where the decays to SM particles also assume 
Set A MSSM input parameters and LO normalisation for the production process; however, 
said contours do not take into account possible SUSY backgrounds. The contour for H" — > bt 
also only takes into account the 36-final state analysis [40]. More detailed studies, including 
46-final states, are ongoing. However, we do not expect major changes in the (m^,tan/3) 
reach for this channel. Similar plots combining the SM and MSSM channels can be drawn 
for the other two MSSM parameter sets. Comparison of the "tb" and "ri/" contours in Fig. 
9 with the analogous discovery regions in [4], which used // = — 200GeV and M2 — 200 GeV 
as inputs, show the former contours to have shrunk somewhat relative to the latter ones, as 
expected since the combined BR(i/^ inos) is larger in relevant parts of the {mA,tan (3) 
plane for Set A inputs than for the inputs of [4]. This shows that the ino decays will reduce 
the rates for the conventional signatures. In particular, relative to a case where the ino 
decay modes are closed (such as when |//|, M2, and sfermion masses are all large) the SM-hke 
discovery regions may be significantly reduced. This makes the search for the 'Si + + f 
signature from H"^ — > inos decays all the more important. 



5 Conclusions 

In summary, we have proven that SUSY decays of charged Higgs bosons can profitably be 
exploited at the LHC in order to detect these important particles. We have done an extensive 
probe of the MSSM parameter space to see where decays of the type xfx^i = li2, 

j = 1,2,3,4) can yield hadronically quiet three lepton (electrons and/or muons) final states. 
Here all tree-level decay chains allowable within the MSSM have been taken into account. 
Coupling such decay chains with top-associated charged Higgs boson production, we selected 
a signature consisting of three hard isolated leptons (electrons and/or muons), three hard 
jets which reconstruct the top quark (with one pair thereof also reconstructing a boson 
and the other bearing a 6-tag) and substantial missing transverse energy. We then performed 
quite realistic MC studies utilising the HERWIG event generator and modelling the CMS 
detector. The hard subprocess used for the signal was gb — > tH~ (and c.c), supplemented by 
initial and final state parton shower and hadronisation, with overall LO normalisation. (All 
backgrounds were generated at the same level of accuracy.) Recent studies [24] have found 
that, in contrast to the negative corrections to said LO production subprocess utilised in the 
past, NLO corrections are in fact positive and may be substantial — depending upon the 
choice of input scales (see [24]) X-factor enhancements of ^ 1.6 may be obtained, comparable 

diminishes as the tan (3 value increases. 

^^The authors caution that this figure is vahd for a specific set of the MSSM inputs M2, and , not in 
general. 
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to or even larger than the corresponding corrections for the leading backgrounds. Inclusion 
of such NLO effects in future signal and background analyses may well expand the discovery 
reach of this channel. 

We found that this '3£ + + signature has the potential to provide coverage over an 
area of the MSSM parameter space roughly corresponding to 250 GcV < m//± < 500 GcV and 
3 < tan/3 < 35. This region covers a substantial portion of parameter space where decays 
into ordinary particles have been shown to be ineffective. However, to this must be added 
the caveat that other MSSM input parameters must be favourable. To wit, a small value for 

and a small to moderate M2 value are essential for having substantial — > xf BRs 
(with M2 > to put more weight on ino decays not including the LSP) and light sleptons 
are crucial for enhancing the leptonic BRs of the inos. Said slepton intermediates may be on- 
or off- mass shell; though of course it is optimal if the two-body on-shell ino decay mode into 
a slepton and a lepton is open, as shown by Fig. 5. Naturally, the actual physical masses of 
the sleptons (selectrons, smuons and the associated sneutrinos) should be less than those of 
an ino pair into which the charged Higgs boson has a significant BR. Depending on the ino 
masses as fixed by the MSSM parameter inputs, this dictates slepton masses of < 160 GeV or 
lower in the discovery regions documented in this work. 

Regions in MSSM parameter space satisfying such criteria tend to be sufficiently close to 
the LEP2 limits and/or to those derived after Run 2 at the Tevatron that such regions should 
be readily accessible to probing by the LHC We have made very few assumptions about 
the underlying SUSY-breaking dynamics associated with some much higher energy scale, and 
hence defined all relevant MSSM input parameters at the EW scale. (The mSUGRA model was 
analysed as a possible GUT benchmark but failed to shown any potential for the considered 
decay channel.) The discovery reach shown in Fig. 9 (for a reasonably favourable choice of 
these parameters) illustrates the possible power of this new channel. 
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Figure 8: 5(j-discovery contours in the tan/? vs. rriA plane for Parameter Sets A, B and C, 
assuming an integrated luminosity of 100 fb~^. The shaded region at the left and bottom of 
each plot is excluded by LEP2 Higgstrahlung {i.e., e^e~ hZ and e^e~ — >■ hA) limits. The 
region in the top plot below the red dotted line is excluded by the LEP2 chargino mass bound. 
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Figure 9: 5cr-discovery contours in the tan/? vs. rriA plane for all charged Higgs channels, both 
SM and MSSM, assuming MSSM inputs as in Parameter Set A and 100 fb~^ of integrated 
luminosity. The area below the red dotted line at the left is excluded by LEP2 Higgstrahlung 
{i.e., e^e~ — > hZ and e+e~ — > hA) limits and the region below the horizontal red dotted line 
is excluded by the LEP2 chargino mass bound. 
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